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Abstract 



We propose in this paper a new model for describing the unification of 
dark energy and dark matter. This new model is a further generalization 
of the generalized Chaplygin gas (GCG) model, thus dubbed new generalized 
Chaplygin gas (NGCG) model. The equation of state of the NGCG is given by 
p = —A(a)/p a , where a is the scale factor and A(a) = — wxAa~ 3 ^ 1+Wx ^ 1+a \ 
We show that the NGCG model is totally dual to an interacting XCDM 
parametrization scenario, in which the interaction between dark energy and 
dark matter is characterized by the constant a. We discuss the cosmological 
consequences led by such an unified dark sectors model. Furthermore, we 
perform a statefinder analysis on this scenario and show the discrimination 
between this scenario and other dark energy models. Finally, a combined 
analysis of the data of Type la supernovae, cosmic microwave background, 
and large scale structure provides a fairly tight constraint on the parameters 
of the NGCG model. 
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I. INTRODUCTION 



Recent observations of Type la supernovae (SNe la) [1] indicate that the expansion of 
the Universe is accelerating at the present time. These results, when combined with the 
observations of cosmic microwave background (CMB) [2] and large scale structure (LSS) [3], 
strongly suggest that the Universe is spatially flat and dominated by an exotic component 
with large negative pressure, referred to as dark energy [4]. The first year result of the 
Wilkinson Microwave Anisotropy Probe (WMAP) shows that dark energy occupies about 
73% of the energy of our Universe, and dark matter about 23%. The usual baryon matter 
which can be described by our known particle theory occupies only about 4% of the total 
energy of the Universe. Although we can affirm that the ultimate fate of the Universe is 
determined by the feature of dark energy, the nature of dark energy as well as its cosmological 
origin remain enigmatic at present. So far the confirmed information about dark energy is 
still limited and can be roughly summarized as the following several items: it is a kind of 
exotic energy form with sufficiently large negative pressure such that drives the Universe 
to undergo a period of accelerating expansion currently; it is spatially homogeneous and 
non-clustering; and it is in small part at the early times while dominates the Universe very 
recently. 

The investigation of the nature of dark energy is an important mission in the modern 
cosmology. Much work has been done on this issue, and there is still a long way to go. 
Currently, the preferred candidates of dark energy are vacuum energy (or cosmological con- 
stant) and dynamical fields. The simplest form of dark energy is the cosmological constant 
A. A tiny positive cosmological constant which can naturally explain the current accelera- 
tion would encounter many theoretical problems, such as the "fine-tuning" problem and the 
"coincidence" problem. Another possible form of dark energy is provided by scalar fields. 
Dark energy can be attributed to the dynamics of a scalar field 0, called quintessence [5], 
which convincingly realize the present accelerated expansion of the Universe by evolving 
slowly down its potential V(<p). The tracker version quintessence is to some extent likely to 
resolve the coincidence problem. It should also be pointed out that the coupled quintessence 
models [6] provide this problem with a more natural solution. However, for quintessence 
models with flat potentials, the quintessence field has to be nearly massless and one thus 
expects radiative corrections to destabilize the ratio between this mass and the other known 
scales of physics. In addition, for the cosmological constant and many quintessence models, 
the event horizon would lead to a potential incompatibility with the string theory. 

Other designs on dark energy roughly include k-essence [7], quiessence (or "X- matter") 
[8], brane world [9], tachyon [10], generalized Chaplygin gas (GCG) [11,12], holographic dark 
energy [13,14], and so forth. The quiessence or X-matter component is simply characterized 
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by a constant, non-positive equation of state wx (wx < — 1/3 is a necessary condition 
to make the Universe accelerate). In general, in a Friedmann- Robertson- Walker (FRW) 
background with the presence of cold dark matter (CDM), an arbitrary but constant wx for 
dark energy from the range (—1, 0) can be achieved by using a scalar field with a hyperbolic 
sine potential [8]. It may be noted that in principle the value of wx may be even less than 
-1. In fact, by fitting the SNe la data in the framework of XCDM (X-matter with CDM), 
the hint for w x < — 1 has been found. Indeed, a study of high-z SNe la [15] finds that 
the equation of state of dark energy has a 99% probability of being < —1 if no priors are 
placed on Q° m . When these SNe results are combined with CMB and 2dFGRS the 95% 
confidence limits on an unevolving equation of state are —1.46 < wx < —0.78 [15,16] which 
is consistent with estimates made by other groups [2,3]. The possibility of wx < —I has 
provoked lots of investigations on the phantom dark energy [17]. The remarkable feature of 
the phantom model is that the Universe will end its life with a "Big Rip" (future singularity) 
within a finite time [18]. On the other hand, we concern here another interesting proposal 
on dark energy, i.e. the dark energy component might be explained by a background fluid 
with an exotic equation of state, the generalized Chaplygin gas model [11]. The striking 
feature of this model is that it allows for a unification of dark energy and dark matter. 
This point can be easily seen from the fact that the GCG behaves as a dust-like matter at 
early times and behaves like a cosmological constant at late stage. This dual role is at the 
heart of the surprising properties of the GCG model. Moreover, the GCG model has been 
successfully confronted with various phenomenological tests involving SNe la data, CMB 
peak locations, gravitational lensing and other observational data [12]. It is remarkable that 
the GCG equation of state has a well defined connection with string and brane theories [19], 
and this gas is the only gas known to admit a supersymmetric generalization [20]. 

In addition, it should be pointed out that the GCG model can be portrayed as a picture 
that cosmological constant type dark energy interacts with cold dark matter. However, since 
the equation of state of dark energy still cannot be determined exactly, the observational data 
show wx is in the range of (—1.46, —0.78), the GCG model should naturally be generalized 
to accommodate any possible X-type dark energy with constant wx- Therefore, we propose 
here a new generalized Chaplygin gas (NGCG) scenario as a scheme for unification of X-type 
dark energy and dark matter. The feature of this new model should also be exhibited in 
that dark sectors are uniformly described by an exotic background fluid and this new gas 
behaves as a dust-like matter at early times and as a X-type dark energy at late times. We 
will show in this paper that this model is a kind of interacting XCDM model, and constrain 
the parameters of this model by using observational data. 

This paper is organized as follows. In Section 2, we introduce the extension version of 
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the generalized Chaplygin gas, namely the NGCG model, to describe the unification of dark 
energy and dark matter, and demonstrate that the NGCG actually is a kind of interacting 
XCDM system. In Section 3, we analyze the NGCG model by means of the statefinder 
parameters. In Section 4, we constrain the parameters of the NGCG model using the SNe 
la, CMB, and LSS data. We give concluding remarks in the final section. 



II. THE NGCG SCENARIO AND THE INTERACTING XCDM 

PARAMETRIZATION 



In this section we introduce the NGCG model. In the framework of FRW cosmology, 
considering an exotic background fluid, the NGCG, described by the equation of state 

A(a) 

Pch = -— , (1) 

where a is a real number and A(a) is a function depends upon the scale factor of the Universe, 
a. We might expect that this exotic background fluid smoothly interpolates between a dust 
dominated phase p ~ a~ 3 and a dark energy dominated phase p ~ a~ 3( - 1+Wx ^ where wx is a 
constant and should be taken as any possible value in the range (—1.46, —0.78). It can be 
expected that the energy density of the NGCG should be elegantly expressed as 



f)ch = Aa S(l+ WX )(l+a) + Sa -3(l+a)l da (2) 



The derivation of the Eq. (2) should be the consequence of substituting the equation of 
state Eq. (1) into the energy conservation equation of the NGCG for an homogeneous and 
isotropic spacetime, this requires the function A(a) to be of the form 

A(a) = -w x Aa- 3{1+Wx ^ 1+ ^ , (3) 

where A is a positive constant, and the other positive constant B appears in Eq. (2) an 
integration constant. One can see explicitly that this model recovers the GCG model as 
the equation-of-state parameter wx taken to be —1, and an ordinary XCDM model can be 
reproduced by taking the parameter a to be zero. The parameter a is called interaction 
parameter of the model as will be shown below. 

The NGCG scenario involves an interacting XCDM picture. For showing this, we first 
decompose the NGCG fluid into two components, one is the dark energy component, and 
the other is the dark matter component, 

PCh = PX + Pdrn ■ (4) 

Note that the pressure of the NGCG fluid is provided only by the dark energy component, 
namely pch = Px- Therefore, the energy density of the dark energy ingredient can be given 
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_ pch _ Aa -3(l+ Wx )(l+ a ) 

PX ~ w x ~ [Aa- 3 ( 1 +^)( 1 + Q ) + Ba-^+*)]Th ' ^ ' 

and then the energy density of the dark matter component can also be obtained 

= Ba-^ 

Pdm [A a -3(n-»»x)(i+«) + 5 a - 3 ( 1 + Q )]^ ' U 

From these expressions one obtains the scaling behavior of the energy densities 

P dm _ ^_ a 3w x {l+a) /j\ 

Px A 

We see explicitly from it that there must exist an energy flow between dark matter and dark 
energy provided that a^O. When a > 0, the transfer direction of the energy flow is from 
dark matter to dark energy; when a < 0, just the reverse. Therefore, it is clear that the 
parameter a characterizes the interaction between dark energy and dark matter. This is the 
reason for that we call a the interaction parameter. 

The parameters A and B can be expressed in terms of current cosmological observables. 
From Eq.(2), it is easy to get 

A + B = pl h0 , (8) 
where i] = 1 + a is used to characterize the interaction for simplicity, thus we have 

A = pl h0 A s , B = pl h0 (l-A s ) , (9) 
where A s is a dimensionless parameter. Using Eqs. (7) and (9), one gets 

a Pxo Q° x 



s 



(10) 



PXO + PdmO 

where the second equality stands for the cosmological model involving the baryon matter 
component. We have assumed here that the space of the Universe is flat. Hence, the NGCG 
energy density can be expressed as 

Pch = Pch a~ 3 [l - Ml - a- 3w ^)]Vv . (11) 

Making use of Eqs. (2), (5), (6), and (11), the energy densities of dark energy and dark 
matter can be re-expressed as 

(12) 



Px = Pxoa- 3(1+ ^ 



0° 



fig 



_3 

Pdm — PdmO a 



x _ _3_ {1 _ a -3^) 



(13) 
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The whole NGCG fluid satisfies the energy conservation, but dark energy and dark 
matter components do not obey the energy conservation separately; they interact with each 
other. We depict this interaction through an energy exchange term Q. The equations of 
motion for dark energy and dark matter can be written as 

p x + 3H(l + w x )p x = Q , (14) 



p dm + 3Hp dm = -Q , (15) 

where dot denotes a derivative with respect to time t, and H = a/a represents the Hubble 
parameter. For convenience we define the effective equations of state for dark energy and 
dark matter through the parameters 

Wx=wx -Wp- X > (16) 

= v§- ■ (17) 

onp dm 

According to the definition of the effective equations of state, the equations of motion for 
dark energy and dark matter can be re-expressed into forms of energy conservation, 

p x + 3H(l + wf)p x = , (18) 

p dm + 3H(l + w d i)p dm = . (19) 

By means of the concrete forms of dark energy and dark matter in NGCG scenario, Eqs. 
(12) and (13), one can obtain 

w (e) _ (y-i)w x (i-n x -ntW w ^ m 
Wx ~ Wx + n x + (i - n x - n°) a ^ ' (2Uj 

( e) = (v - i)w x n° x 

dm n° x + (i - n° x - n° b )a 3w ^ ' 1 ; 

Now we switch to discuss the cosmological evolution. Consider a spatially flat FRW 
Universe with baryon matter component pb and NGCG fluid pch, the Friedmann equation 
reads 

3M 2 P H 2 = p Ch + p b , (22) 
where Mp is the reduced Planck mass. The Friedmann equation can also be expressed as 
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H(a) = H E(a) , 



(23) 



where 



E(a) 



1 - 



n 



X 



i-ng 



(1 - a 



-3w x r)\ 



1/r, 



1/2 



+ n h a~ 3 



(24) 



Then, the fractional energy densities of various components can be easily obtained 



tt x = n x E- 2 a- 3( - 1+Wxv) 



(25) 



Q dm = (i - n° x - n° b )E- 2 a 



2„-3 



T-a 



-3u>x7A 



(26) 



(27) 



So far we see clearly that the NGCG model is totally dual to a coupled dark energy 
scenario [21,22], namely an interacting XCDM parametrization. It is remarkable that the 
interaction between dark energy and dark matter can be interpreted as arising from the 
time variation of the mass of dark matter particles. The GCG model is a special case in the 
NGCG model corresponding to w x — — 1, thus the GCG model is actually an interacting 
ACDM model. Fig.l and Fig.2 illustrate examples of the density evolution in the NGCG 
model. The current density parameters used in the plots are VP dm = 0.25, Vt x = 0.7, and 
f2° = 0.05. In Fig.l, we show the cases having the common equation-of-state parameter 
wx = —1-2, while the interaction parameters a are taken to be 0, 0.5, and 0.8, respectively. 
Note that the a = case corresponds to a normal phantom model with constant wx- In this 
example we see the role the interaction parameter a plays in the model. The transfer energy 
flows from dark matter to dark energy when a > 0; the larger a leads to the stronger energy 
flow; density of baryon component is also affected evidently by the interaction between dark 
energy and dark matter. In Fig.2, we depict the cases with common interaction parameter 
a = 0.5, and the equation-of-state parameters w x are taken to be —1, —0.8, and —1.2, 
respectively. Here w x = — 1 case corresponds exactly to the GCG model. The effect of the 
parameter wx in the NGCG scenario is also evident as we see in this example. 
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FIG. 1. The evolution of the density parameters for various components Qx, ^dmi an d 
Note that Qch = ^x + ^dm- The current density parameters used in the plot are £l dm = 0.25, 
Q x = 0.7, and = 0.05. In this case, we fix wx and vary a. 




2 3 

redshift z 

FIG. 2. The evolution of the density parameters for various components tlx, ^dm, and 
Note that f^ch = ^x + ^dm- The current density parameters used in the plot are fi[j m = 0.25, 
Q x = 0.7, and = 0.05. In this case, we fix a and vary wx- 



Let us now discuss the cosmological consequences led by the NGCG model and compare 
the cosmological quantities in the NGCG cosmology with those of some special cases such as 
ACDM and GCG. Firstly, we regard the Hubble parameter H which evaluates the expansion 
rate of the Universe. In Fig. 3 we plot the Hubble parameter of the NGCG model in units of 
H\cdm as a function of redshift z range from to 5. The current density parameters used in 
the plot of Fig. 3 are the same as used in Figs.l and 2. The model parameters are divided into 
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two groups, a = and a = 0.5, both including wx = — 1, —0.8, and —1.2. It can be seen 
from Fig.3 that the NGCG model degenerates to the XCDM when the parameter a takes 0; 
the cases of wx > — 1 and wx < — 1 make H larger than and less than H\ C dm, respectively, 
during the cosmological evolution. The introducing of the interaction parameter a makes H 
be larger than i?ACDM evidently at early times; while it is interesting to see that the value 
of H/H AC dm can cross 1 in the case of w x < — 1 in recent period. The acceleration of the 
Universe is evaluated by the deceleration parameter q = —d/aH 2 . Omitting the radiation 
component, the deceleration parameter can be expressed as 

1=\ + l^x , (28) 

where Qx is given by (25). The evolution of the deceleration parameter q is depicted in 
Fig.3 for selected parameter sets. The current density parameters are taken to be the same 
as above figures. The influence coming from the interaction a and equation of state of dark 
energy wx can be seen clearly in this figure. We notice that a positive a makes the redshift 
of acceleration/deceleration transition (q(z T ) = 0) shift to a smaller value; while the values 
of zt are nearly degenerate under the same a as shown in this example. 
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FIG. 3. The evolution of the Hubble parameter in units of -£/acdm(^). The current density 
parameters are taken to be Q® m = 0.25, Q x = 0.7, and = 0.05. 
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a=0,W x =-1 (ACDM) 
a=0.5,w x =-1 (GCG) 
a=0.5,w x =-0.8 
a=0.5,w =-1.2 
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FIG. 4. The evolution of the deceleration parameter q(z). The current density parameters are 
taken to be tt° dm = 0.25, Q° x = 0.7, and fig = 0.05. 



III. STATEFINDER DIAGNOSTIC 

Since more and more dark energy models have been constructed for interpreting or 
describing the cosmic acceleration, the problem of discriminating between the various con- 
tenders becomes very important. In order to be able to differentiate between those competing 
cosmological scenarios involving dark energy, a sensitive and robust diagnostic for dark en- 
ergy models is a must. For this purpose a diagnostic proposal that makes use of parameter 
pair {r, s}, the so-called "statefinder" , was introduced by Sahni et al. [23]. The statefinder 
probes the expansion dynamics of the Universe through higher derivatives of the scale factor 
a and is a natural companion to the deceleration parameter q which depends upon a. The 
statefinder pair {r, s} is defined as follows 

a r — 1 . . 

(29) 



a# 3 ' 3(q- 1/2) ' 

The statefinder is a "geometrical" diagnostic in the sense that it depends upon the scale 
factor and hence upon the metric describing space-time. 

Trajectories in the s — r plane corresponding to different cosmological models exhibit 
qualitatively different behaviors. The spatially flat ACDM scenario corresponds to a fixed 
point in the diagram 

{s,r} ={0,1}. (30) 

ACDM 

Departure of a given dark energy model from this fixed point provides a good way of estab- 
lishing the "distance" of this model from ACDM [23,24]. As demonstrated in Refs. [23-29] 
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the statefinder can successfully differentiate between a wide variety of dark energy models 
including the cosmo logical constant, quintessence, quintom, the Chaplygin gas, braneworld 
models, holographic dark energy and interacting dark energy models. We can clearly iden- 
tify the "distance" from a given dark energy model to the ACDM scenario by using the r(s) 
evolution diagram. 

The current location of the parameters s and r in these diagrams can be calculated 
in models, and on the other hand it can also be extracted from data coming from SNAP 
(SuperNovae Acceleration Probe) type experiments [23,24]. Therefore, the statefinder di- 
agnostic combined with future SNAP observations may possibly be used to discriminate 
between different dark energy models. For example, as shown in Ref. [24], by carrying out 
a maximum likelihood analysis which combines the statefinder diagnostic with realistic ex- 
pectations from the SNAP experiment, the averaged-over-redshift statefinder pair {s,r} is 
convincingly demonstrated to be useful diagnostic tool in successfully differentiating between 
the cosmological constant and dynamical models of dark energy. In this section we apply 
the statefinder diagnostic to the NGCG model. 

In what follows we will calculate the rangefinder parameters for the NGCG model 
and plot the evolution trajectories of the model in the statefinder parameter-plane. The 
statefinder parameters can be expressed in terms of the total energy density p and the total 
pressure p in the Universe: 

r = l + 9(p + p) g, S = (P±AP. (31) 
2p p p p 

The total energy of the Universe is conserved, so we have p = —3H(p + p). Since the dust 
matter does not have pressure, the total pressure of the cosmic fluids is provided only by 
dark energy component, p = px = WxPx- Then making use of p — —3H(p + p) and 
px = —3H(1 + w^x)px, we can get the concrete expression of the statefinder parameters 

r = l + ^w x n x (l + w% ) ) , s = l+w^. (32) 

Here w^ x and Q x are given by (20) and (25), respectively. Though the relationship between 
statefinder parameters r and s, namely the function r(s), might be derived analytically in 
principle, we do not give the expression here due to the complexity of the formula. Making 
the redshift z — 1/a — 1 vary in an enough large range involving far future and far past, e.g. 
from —1 to 5, one can easily get the evolution trajectories in the statefinder s — r plane of 
this model. Selected curves of r(s) are plotted in Fig.5 and Fig. 6. In Fig. 5, we fix a = 0.5 
and vary wx as —0.8, —1, and —1.2, respectively. In Fig.6, we fix wx = —1.2 and vary a as 
0, ±0.2, ±0.5, and ±0.8, respectively. Other parameters are taken as the same as previous 
figures. In these two figures, dots locate the today's values of the statefinder parameters 
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(so,r ) and arrows denote the evolution directions of the statefmder trajectories r(s). The 
ACDM model locates at (0, 1) in the s — r plane also denoted as a dot. 




-0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 



S 

FIG. 5. The statefmder r(s) evolution diagram. Dots locate the today's values of the statefmder 
parameters (so> r o) an d arrows denote the evolution directions of the statefmder trajectories r(s). 
The ACDM model locates at the fixed point (0, 1). In this case, we fix a = 0.5 and vary wx as 
—0.8, —1, and —1.2, respectively. 



i 1 . 1 . 1 1 1 1 1 1 1 . 1 . 1 r 




-1.4 -1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 

S 

FIG. 6. The statefmder r(s) evolution diagram. Dots locate the today's values of the statefmder 
parameters (so,?"o) an d arrows denote the evolution directions of the statefmder trajectories r(s). 
The ACDM model locates at the fixed point (0, 1). In this case, we fix wx = —1.2 and vary a as 
0, ±0.2, ±0.5, and ±0.8, respectively. 
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The statefinder diagnostic can discriminate between various dark energy models effec- 
tively. Different cosmological models involving dark energy exhibit qualitatively different 
evolution trajectories in the s — r plane. For example, the ACDM scenario corresponds 
to the fixed point s = 0, r = 1 as shown in (30), and the SCDM (standard cold dark 
matter) scenario corresponds to the point s — 1, r = 1. For the "quiessence" (XCDM) 
models, the trajectories are some vertical segments, i.e. r decreases monotonically from 1 
to 1 + |u>x(l + w x ) while s remains constant at 1 +w x [23,24]. The quintessence (inverse 
power law) tracker models have typical trajectories similar to arcs of an upward parabola 
lying in the regions s > 0, r < 1 [23,24]. The holographic dark energy scenario (c = 1 
case), as shown in [27], commences its evolution from s — 2/3, r — 1, through an arc 
segment, and ends it at the ACDM fixed point (s — 0, r = 1) in the future. The coupled 
quintessence models and Quintom models exhibit more complicated trajectories as shown in 
Refs. [25,29]. Now from the Figs. 5 and 6 of this paper, we can see the statefinder trajecto- 
ries of the NGCG model. In Fig. 5 we see the cases under the fixed a, where the GCG model 
(wx = — 1) exhibits a complete downward parabola, while the general cases (w x 7^—1) cor- 
respond to some broken parabolas. The statefinder trajectory commences its evolution from 
s = 1 + (1 + a)wx, r — lait— > to s = 1 + wx, r = 1 + §u>x(l + wx) at t — > oo, through 
the curve. Today's statefinder point locates at s = 1 + w xo-> r o = 1 + ^ w x^°x(^ + w xo)-> 
where w%l = w x + aw x (l - tt° x - - fig). The "distance" from the NGCG model 

to the ACDM scenario can be measured directly in the statefinder plane. Note that under 
a positive a, the cases of w x < —1 never arrive at the ACDM fixed point; the GCG case 
(wx = — 1) ends at the ACDM fixed point; while the cases of wx > — 1 have passed through 
this fixed point. Fig.6 displays the cases under a fixed wx- We show here wx = —1-2, a 
phantom. Trajectories correspond to zero, positive, as well as negative values of a are all 
displayed in this diagram to depict a complete statefinder diagnostic. It is interesting to see 
that the trajectories can pass through the ACDM fixed point under a phantom case when 
a < 0. This is because a negative a makes dark energy component transfer energy flow 
to dark matter component. We notice that the normal phantom case (a = 0) evolves its 
trajectory along a vertical segment. Comparing with the quiessence case [23,24], the phan- 
tom case reposes on the left of ACDM point, namely the region s < 0, r > 1, and evolves 
upwards; while the quiessence case reposes on the right of the ACDM point, namely the 
region s > 0, r < 1, and evolves downwards. Interestingly, under a fixed wx, the present 
statefinder points as well as the ACDM fixed point locate on a straight line. This is because 
when wx is fixed, the relationship between r and s is linear, r = 1 + ^wx^x s o- 
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IV. OBSERVATIONAL CONSTRAINTS FROM SNE IA, CMB, AND LSS DATA 



In this section we will derive the constraints on the NGCG model from current available 
observational data. It should be mentioned that the interacting XCDM parametrization 
scenario has been tested by the recent Type la supernovae data [21]. The results show 
that the SNe la data favor a negative coupling and an equation of state wx < — 1, namely 
a negatively coupled phantom dark energy. However, as we know, the supernovae data 
alone are not sufficient to constrain dark energy models strictly (see e.g. the analysis in 
Ref. [14]). Therefore, to obtain more tight constraints on dark energy models, one should 
need additional data provided by other astronomical observations to be necessary and useful 
complements to the SNe data. It has been demonstrated that some observational quantities 
irrelevant to H are very suitable to play this role [30]. Such quantities and data can be 
found in the probes of CMB and LSS [14,30-33]. In what follows we perform a combined 
analysis of SNe la, CMB, and LSS on the constraints of the NGCG model. We use a x 2 
statistic 



where x!n> Xcmb an d Xlss are contributions from SNe la, CMB, and LSS data, respectively. 
It is well known that the acceleration of the Universe is found by the Type la supernovae 
observations, where the concept of the luminosity distance plays a very important role. The 
luminosity distance of a light source is defined in such a way as to generalize to an expanding 
and curved space the inverse-square law of brightness valid in a static Euclidean space, 



where C is the absolute luminosity which is a known value for the standard candle SNe la, T 
is the measured flux. The Hubble distance cHq 1 = 2997. 9/iT 1 Mpc. The Type la supernova 
observations directly measure the apparent magnitude m of a supernova and its red-shift 
z. The apparent magnitude m is related to the luminosity distance di of the supernova 
through 



where M is the absolute magnitude which is believed to be constant for all Type la super- 
novae. In our analysis, we take the 157 gold data points listed in Riess et al. [16] which 
includes recent new 14 high redshift SNe (gold) data from the HST/GOODS program. The 
X 2 function determined by SNe la observations is 



Z Z I Z I z 

X — XSN XCMB XlSS > 



(33) 




(34) 



m(z) = M + 5 log 10 (rf L (^)/Mpc) + 25 , 



(35) 



14 



157 




(36) 



i=i 



where the extinction-corrected distance moduli /i(z) is defined as /i(z) = m(z) — M, and 
<7j is the total uncertainty in the observation. Following the Ref. [21], we fix = 0.05 
in the computation for simplification. Hence, the computation is carried out in a four- 
dimensional space, for the four parameters P = (r),wx, h,Q dm ). For the CMB, we use only 
the measurement of the CMB shift parameter [34], 



where = £l dm + and z dcc = 1089 [2]. Note that this quantity is irrelevant to the 
parameter H such that provides robust constraint on the dark energy model. The results 
from CMB data correspond to TZ = 1.716 ±0.062 (given by WMAP, CBI, ACBAR) [2,35]. 
We include the CMB data in our analysis by adding Xcmb = [(^ ~~ T^o)/c-r] 2 (see e.g. Refs. 
[30-32]), where 1Z is computed by the NGCG model using equation (37). The only large 
scale structure information we use is the parameter A measured by SDSS [36], defined by 



where z\ = 0.35. Also, we find that this quantity is independent of Hq either, thus can 
provide another robust constraint on the model. The SDSS gives the measurement data 
[36] Ao = 0.469 ± 0.017. We also include the LSS constraint in our analysis by adding 
Xlss = [(^ — Ao)/<Ja} 2 (see e.g. Refs. [14,33]), where A is computed by the NGCG model 
using equation (38). Note that we have chosen to use only the most conservative and robust 
information, TZ and A, from CMB and LSS observations. These measurements we use do 
not depend on the Hubble constant Hq, thus are useful complements to the SNe data. It 
is remarkable that the likelihood analysis scheme we employ here is very economical and 
efficient due to that it does not make use of all the information available in CMB and LSS 
but can provide fairly good constraints on dark energy models [14,33]. 

We now analyze the probability distribution of rj and wx in the NGCG model. The like- 
lihood of these two parameters is determined by minimizing over the "nuisance" parameters 



where the integral is over a large enough range of h and Vt dm to include almost all the 
probability. We now compute C(i], wx) on a two-dimensional grid spanned by r] and wx- The 
68.3%, 95.4%, and 99.7% (namely 1, 2, and 3 a) confidence contours consist of points where 
the likelihood equals e~ 2 - 31 / 2 , e ~ 6 - 18 / 2 , and e - 1L83 / 2 Q f the maximum value of the likelihood, 




(37) 




(38) 




(39) 
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respectively. Fig. 7 shows our main results, the contours of la, 2a, and 3a confidence levels 
in the w x — V plane. The 1 a fit values for the model parameters are: w x = — 0-98tg;2o 
and r] = 1.061&16, and the minimum value of x 2 in the four dimensional parameter space is: 
Xmin = 167.29. We see clearly that the combined analysis of SNe la, CMB, and LSS data 
provides a fairly tight constraint on the NGCG model. It is remarkable that the best fit 
happens at the vicinity of the cosmological constant, even though w x is slightly larger than 
— 1 and i] is mildly larger than 1 (i.e. a slightly larger than 0). This means that within the 
framework of the NGCG model the real form of dark energy at the maximum probability 
is the near cosmological constant according to the joint analysis of SNe+CMB+LSS data. 
However, the analysis results still accommodate the existence likelihood of the "X-matter" 
and the interaction between dark energy and dark matter. In 1 a range, wx G (—1.18, —0.83) 
and a G (—0.1,0.26). This implies that the probabilities of that dark energy behaves as 
quintessence-like form and phantom-like form are roughly equal, and the probabilities that 
the energy flow streams from dark energy to dark matter and the reverse are also roughly 
equal. One-dimensional likelihood distribution functions for wx and rj are shown in Fig. 8 and 
Fig. 9, respectively. It is very clear that the original Chaplygin gas model, a — 1 (or rj = 2) 
and wx = — 1, is totally ruled out by the observational data at 99.7% confidence level. In 
addition, it should be pointed out that when we fix rj = 2 and let w x free, the NGCG model 
will be identified as the so-called variable Chaplygin gas (VCG) model proposed in Ref. [37] 
(see also Ref. [38]). The joint analysis of SNe+CMB+LSS also rules out this probability. It 
is hopeful that the future precise data will provide more strong evidences to judge whether 
the dark energy is the cosmological constant and whether dark energy and dark matter are 
in unification. 
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FIG. 7. Confidence level contours of 68.3%, 95.4% and 99.7% in the (r),w x ) plane. The 1 a fit 
values for the model parameters are: wx = — 0.98lo2o an d V = l-06lg i6> an d ^ ne minimum value 
of x 2 hi the f° ur dimensional parameter space is: Xmm = 167.29. 
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FIG. 8. One-dimensional probability distribution for wx- 
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FIG. 9. One-dimensional probability distribution for 77. 



V. CONCLUDING REMARKS 

The Chaplygin gas model is a proposal to describe dark energy and dark matter as a 
unified fluid, the Chaplygin gas, characterized by an exotic equation of state p = —A/p, 
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where A is a positive constant. Since this original Chaplygin gas has been ruled out by the 
observations, a generalization of the Chaplygin equation of state p = —A/p a was considered, 
by introducing a free parameter a. The generalized Chaplygin gas model is also regarded 
as a unification of dark energy and dark matter. The reason is that the GCG behaves as a 
dust-like matter at early stage and as a cosmological constant at late stage. That is to say 
that the GCG model admits that the Universe will be dominated by a cosmological constant 
and thus enter into a de Sitter phase in the future. However, we can not hitherto affirm 
whether the dark energy is a tiny positive cosmological constant. Therefore, the scheme for 
unification of dark energy and dark matter should accommodate other forms of dark energy 
such as quintessence-like and phantom-like dark energy. This is the motivation for us to 
further generalize the GCG model. 

We propose in this paper a new model as a scheme for the unification of dark energy 
and dark matter. This new model is a further extension version of the GCG model, thus 
dubbed new generalized Chaplygin gas model. This further generalization is implemented 
by introducing another free parameter wx to make the constant A in the GCG equation of 
state become a scale factor dependent function A(a). In order to implement the interpolation 
between a dust dominated Universe and an X-matter dominated Universe, the unique choice 
of A(a) is A(a) = — w x Aa^ 1+Wx ^ 1+a \ Through a two-fluid decomposition, we show that 
the NGCG model is totally equivalent to an interacting XCDM parametrization scenario, in 
which the interaction between dark energy and dark matter is characterized by the constant 
a. We discuss the cosmological consequences led by such an unified dark sectors model. 
Furthermore, a statefinder diagnostic is performed on this scenario and the discrimination 
between this scenario and other dark energy models is shown. Finally, a combined analysis of 
the data of SNe la, CMB, and LSS is used to constrain the parameters of the NGCG model. 
The fit result shows that the joint analysis can provide a considerably tight constraint on 
the NGCG model. According to the observational test, the best fit happens at the vicinity 
of the ACDM. We hope that the future precise data will provide more strong evidences to 
judge whether the dark energy is the cosmological constant and whether dark energy and 
dark matter can be unified into one component. Also, it would be interesting to investigate 
the evolution of density perturbations and the structure formation in the NGCG scenario. 
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